Testing for differential abundance in mass cytometry data.

When comparing biological conditions using mass cytometry data, a key challenge is to identify cellular populations that change in abundance. Here, we present a computational strategy for detecting 'differentially abundant' populations by assigning cells to hyperspheres, testing for significant differences between conditions and controlling the spatial false discovery rate. Our method (http://bioconductor.org/packages/cydar) outperforms other approaches in simulations and finds novel patterns of differential abundance in real data.This work was supported by Cancer Research UK (core funding to J.C.M., award no. A17197), the University of Cambridge and Hutchison Whampoa Limited. J.C.M. was also supported by core funding from EMBL

Exploring movement patterns and changing distributions of baleen whales in the western North Atlantic using a decade of passive acoustic data

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Davis, G. E., Baumgartner, M. F., Corkeron, P. J., Bell, J., Berchok, C., Bonnell, J. M., Thornton, J. B., Brault, S., Buchanan, G. A., Cholewiak, D. M., Clark, C. W., Delarue, J., Hatch, L. T., Klinck, H., Kraus, S. D., Martin, B., Mellinger, D. K., Moors-Murphy, H., Nieukirk, S., Nowacek, D. P., Parks, S. E., Parry, D., Pegg, N., Read, A. J., Rice, A. N., Risch, D., Scott, A., Soldevilla, M. S., Stafford, K. M., Stanistreet, J. E., Summers, E., Todd, S., & Van Parijs, S. M. Exploring movement patterns and changing distributions of baleen whales in the western North Atlantic using a decade of passive acoustic data. Global Change Biology, (2020): 1-30, doi:10.1111/gcb.15191.Six baleen whale species are found in the temperate western North Atlantic Ocean, with limited information existing on the distribution and movement patterns for most. There is mounting evidence of distributional shifts in many species, including marine mammals, likely because of climate‐driven changes in ocean temperature and circulation. Previous acoustic studies examined the occurrence of minke (Balaenoptera acutorostrata ) and North Atlantic right whales (NARW; Eubalaena glacialis ). This study assesses the acoustic presence of humpback (Megaptera novaeangliae ), sei (B. borealis ), fin (B. physalus ), and blue whales (B. musculus ) over a decade, based on daily detections of their vocalizations. Data collected from 2004 to 2014 on 281 bottom‐mounted recorders, totaling 35,033 days, were processed using automated detection software and screened for each species' presence. A published study on NARW acoustics revealed significant changes in occurrence patterns between the periods of 2004–2010 and 2011–2014; therefore, these same time periods were examined here. All four species were present from the Southeast United States to Greenland; humpback whales were also present in the Caribbean. All species occurred throughout all regions in the winter, suggesting that baleen whales are widely distributed during these months. Each of the species showed significant changes in acoustic occurrence after 2010. Similar to NARWs, sei whales had higher acoustic occurrence in mid‐Atlantic regions after 2010. Fin, blue, and sei whales were more frequently detected in the northern latitudes of the study area after 2010. Despite this general northward shift, all four species were detected less on the Scotian Shelf area after 2010, matching documented shifts in prey availability in this region. A decade of acoustic observations have shown important distributional changes over the range of baleen whales, mirroring known climatic shifts and identifying new habitats that will require further protection from anthropogenic threats like fixed fishing gear, shipping, and noise pollution.We thank Chris Pelkie, David Wiley, Michael Thompson, Chris Tessaglia‐Hymes, Eric Matzen, Chris Tremblay, Lance Garrison, Anurag Kumar, John Hildebrand, Lynne Hodge, Russell Charif, Kathleen Dudzinski, and Ann Warde for help with project planning, field work support, and data management. For all the support and advice, thanks to the NEFSC Protected Species Branch, especially the passive acoustics group, Josh Hatch, and Leah Crowe. We thank the field and crew teams on all the ships that helped in the numerous deployments and recoveries. This research was funded and supported by many organizations, specified by projects as follows: data recordings from region 1 were provided by K. Stafford (funding: National Science Foundation #NSF‐ARC 0532611). Region 2 data: D. K. Mellinger and S. Nieukirk, National Oceanic and Atmospheric Administration (NOAA) PMEL contribution #5055 (funding: NOAA and the Office of Naval Research #N00014–03–1–0099, NOAA #NA06OAR4600100, US Navy #N00244‐08‐1‐0029, N00244‐09‐1‐0079, and N00244‐10‐1‐0047). Region 3A data: D. Risch (funding: NOAA and Navy N45 programs). Region 3 data: H. Moors‐Murphy and Fisheries and Oceans Canada (2005–2014 data), and the Whitehead Lab of Dalhousie University (eastern Scotian Shelf data; logistical support by A. Cogswell, J. Bartholette, A. Hartling, and vessel CCGS Hudson crew). Emerald Basin and Roseway Basin Guardbuoy data, deployment, and funding: Akoostix Inc. Region 3 Emerald Bank and Roseway Basin 2004 data: D. K. Mellinger and S. Nieukirk, NOAA PMEL contribution #5055 (funding: NOAA). Region 4 data: S. Parks (funding: NOAA and Cornell University) and E. Summers, S. Todd, J. Bort Thornton, A. N. Rice, and C. W. Clark (funding: Maine Department of Marine Resources, NOAA #NA09NMF4520418, and #NA10NMF4520291). Region 5 data: S. M. Van Parijs, D. Cholewiak, L. Hatch, C. W. Clark, D. Risch, and D. Wiley (funding: National Oceanic Partnership Program (NOPP), NOAA, and Navy N45). Region 6 data: S. M. Van Parijs and D. Cholewiak (funding: Navy N45 and Bureau of Ocean and Energy Management (BOEM) Atlantic Marine Assessment Program for Protected Species [AMAPPS] program). Region 7 data: A. N. Rice, H. Klinck, A. Warde, B. Martin, J. Delarue, and S. Kraus (funding: New York State Department of Environmental Conservation, Massachusetts Clean Energy Center, and BOEM). Region 8 data: G. Buchanan, and K. Dudzinski (funding: New Jersey Department of Environmental Protection and the New Jersey Clean Energy Fund) and A. N. Rice, C. W. Clark, and H. Klinck (funding: Center for Conservation Bioacoustics at Cornell University and BOEM). Region 9 data: J. E. Stanistreet, J. Bell, D. P. Nowacek, A. J. Read, and S. M. Van Parijs (funding: NOAA and US Fleet Forces Command). Region 10 data: L. Garrison, M. Soldevilla, C. W. Clark, R. A. Chariff, A. N. Rice, H. Klinck, J. Bell, D. P. Nowacek, A. J. Read, J. Hildebrand, A. Kumar, L. Hodge, and J. E. Stanistreet (funding: US Fleet Forces Command, BOEM, NOAA, and NOPP). Region 11 data: C. Berchok as part of a collaborative project led by the Fundacion Dominicana de Estudios Marinos, Inc. (Dr. Idelisa Bonnelly de Calventi; funding: The Nature Conservancy [Elianny Dominguez]) and D. Risch (funding: World Wildlife Fund, NOAA, and Dutch Ministry of Economic Affairs)

Building a model: developing genomic resources for common milkweed (Asclepias syriaca) with low coverage genome sequencing

<p>Abstract</p> <p>Background</p> <p>Milkweeds (<it>Asclepias </it>L.) have been extensively investigated in diverse areas of evolutionary biology and ecology; however, there are few genetic resources available to facilitate and compliment these studies. This study explored how low coverage genome sequencing of the common milkweed (<it>Asclepias syriaca </it>L.) could be useful in characterizing the genome of a plant without prior genomic information and for development of genomic resources as a step toward further developing <it>A. syriaca </it>as a model in ecology and evolution.</p> <p>Results</p> <p>A 0.5× genome of <it>A. syriaca </it>was produced using Illumina sequencing. A virtually complete chloroplast genome of 158,598 bp was assembled, revealing few repeats and loss of three genes: <it>accD, clpP</it>, and <it>ycf1</it>. A nearly complete rDNA cistron (18S-5.8S-26S; 7,541 bp) and 5S rDNA (120 bp) sequence were obtained. Assessment of polymorphism revealed that the rDNA cistron and 5S rDNA had 0.3% and 26.7% polymorphic sites, respectively. A partial mitochondrial genome sequence (130,764 bp), with identical gene content to tobacco, was also assembled. An initial characterization of repeat content indicated that Ty1/<it>copia</it>-like retroelements are the most common repeat type in the milkweed genome. At least one <it>A. syriaca </it>microread hit 88% of <it>Catharanthus roseus </it>(Apocynaceae) unigenes (median coverage of 0.29×) and 66% of single copy orthologs (COSII) in asterids (median coverage of 0.14×). From this partial characterization of the <it>A. syriaca </it>genome, markers for population genetics (microsatellites) and phylogenetics (low-copy nuclear genes) studies were developed.</p> <p>Conclusions</p> <p>The results highlight the promise of next generation sequencing for development of genomic resources for any organism. Low coverage genome sequencing allows characterization of the high copy fraction of the genome and exploration of the low copy fraction of the genome, which facilitate the development of molecular tools for further study of a target species and its relatives. This study represents a first step in the development of a community resource for further study of plant-insect co-evolution, anti-herbivore defense, floral developmental genetics, reproductive biology, chemical evolution, population genetics, and comparative genomics using milkweeds, and <it>A. syriaca </it>in particular, as ecological and evolutionary models.</p

Lewis Base Mediated β-Elimination and Lewis Acid Mediated Insertion Reactions of Disilazido Zirconium Compounds

The reactivity of a series of disilazido zirconocene complexes is dominated by the migration of anionic groups (hydrogen, alkyl, halide, OTf) between the zirconium and silicon centers. The direction of these migrations is controlled by the addition of two-electron donors (Lewis bases) or two-electron acceptors (Lewis acids). The cationic nonclassical [Cp2ZrN(SiHMe2)2]+ ([2]+) is prepared from Cp2Zr{N(SiHMe2)2}H (1) and B(C6F5)3 or [Ph3C][B(C6F5)4], while reactions of B(C6F5)3 and Cp2Zr{N(SiHMe2)2}R (R = Me (3), Et (5), n-C3H7 (7), CH═CHSiMe3 (9)) provide a mixture of [2]+ and [Cp2ZrN(SiHMe2)(SiRMe2)]+. The latter products are formed through B(C6F5)3 abstraction of a β-H and R group migration from Zr to the β-Si center. Related β-hydrogen abstraction and X group migration reactions are observed for Cp2Zr{N(SiHMe2)2}X (X = OTf (11), Cl (13), OMe (15), O-i-C3H7 (16)). Alternatively, addition of DMAP (DMAP = 4-(dimethylamino)pyridine) to [2]+ results in coordination to a Si center and hydrogen migration to zirconium, giving the cationic complex [Cp2Zr{N(SiHMe2)(SiMe2DMAP)}H]+ ([19]+). Related hydrogen migration occurs from [Cp2ZrN(SiHMe2)(SiMe2OCHMe2)]+ ([18]+) to give [Cp2Zr{N(SiMe2DMAP)(SiMe2OCHMe2)}H]+ ([22]+), whereas X group migration is observed upon addition of DMAP to [Cp2ZrN(SiHMe2)(SiMe2X)]+ (X = OTf ([12]+), Cl ([14]+)) to give [Cp2Zr{N(SiHMe2)(SiMe2DMAP)}X]+ (X = OTf ([26]+), Cl ([20]+)). The species involved in these transformations are described by resonance structures that suggest β-elimination. Notably, such pathways are previously unknown in early metal amide chemistry. Finally, these migrations facilitate direct Si–H addition to carbonyls, which is proposed to occur through a pathway that previously had been reserved for later transition metal compounds

T cell cytolytic capacity is independent of initial stimulation strength.

How cells respond to myriad stimuli with finite signaling machinery is central to immunology. In naive T cells, the inherent effect of ligand strength on activation pathways and endpoints has remained controversial, confounded by environmental fluctuations and intercellular variability within populations. Here we studied how ligand potency affected the activation of CD8+ T cells in vitro, through the use of genome-wide RNA, multi-dimensional protein and functional measurements in single cells. Our data revealed that strong ligands drove more efficient and uniform activation than did weak ligands, but all activated cells were fully cytolytic. Notably, activation followed the same transcriptional pathways regardless of ligand potency. Thus, stimulation strength did not intrinsically dictate the T cell-activation route or phenotype; instead, it controlled how rapidly and simultaneously the cells initiated activation, allowing limited machinery to elicit wide-ranging responses

The First Post-Kepler Brightness Dips of KIC 8462852

We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in October 2015, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1-2.5% dips, named "Elsie," "Celeste," "Skara Brae," and "Angkor", which persist on timescales from several days to weeks. Our main results so far are: (i) there are no apparent changes of the stellar spectrum or polarization during the dips; (ii) the multiband photometry of the dips shows differential reddening favoring non-grey extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale <<1um, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term "secular" dimming, which may be caused by independent processes, or probe different regimes of a single process

Many Labs 5:Testing pre-data collection peer review as an intervention to increase replicability

Replication studies in psychological science sometimes fail to reproduce prior findings. If these studies use methods that are unfaithful to the original study or ineffective in eliciting the phenomenon of interest, then a failure to replicate may be a failure of the protocol rather than a challenge to the original finding. Formal pre-data-collection peer review by experts may address shortcomings and increase replicability rates. We selected 10 replication studies from the Reproducibility Project: Psychology (RP:P; Open Science Collaboration, 2015) for which the original authors had expressed concerns about the replication designs before data collection; only one of these studies had yielded a statistically significant effect (p < .05). Commenters suggested that lack of adherence to expert review and low-powered tests were the reasons that most of these RP:P studies failed to replicate the original effects. We revised the replication protocols and received formal peer review prior to conducting new replication studies. We administered the RP:P and revised protocols in multiple laboratories (median number of laboratories per original study = 6.5, range = 3?9; median total sample = 1,279.5, range = 276?3,512) for high-powered tests of each original finding with both protocols. Overall, following the preregistered analysis plan, we found that the revised protocols produced effect sizes similar to those of the RP:P protocols (?r = .002 or .014, depending on analytic approach). The median effect size for the revised protocols (r = .05) was similar to that of the RP:P protocols (r = .04) and the original RP:P replications (r = .11), and smaller than that of the original studies (r = .37). Analysis of the cumulative evidence across the original studies and the corresponding three replication attempts provided very precise estimates of the 10 tested effects and indicated that their effect sizes (median r = .07, range = .00?.15) were 78% smaller, on average, than the original effect sizes (median r = .37, range = .19?.50)

In vitro and in vivo functions of SARS-CoV-2 infection-enhancing and neutralizing antibodies

SARS-CoV-2 neutralizing antibodies (NAbs) protect against COVID-19. A concern regarding SARS-CoV-2 antibodies is whether they mediate disease enhancement. Here, we isolated NAbs against the receptor-binding domain (RBD) and the N-terminal domain (NTD) of SARS-CoV-2 spike from individuals with acute or convalescent SARS-CoV-2 or a history of SARS-CoV infection. Cryo-electron microscopy of RBD and NTD antibodies demonstrated function-specific modes of binding. Select RBD NAbs also demonstrated Fc receptor-g (FcgR)-mediated enhancement of virus infection in vitro, while five non-neutralizing NTD antibodies mediated FcgR-independent in vitro infection enhancement. However, both types of infection-enhancing antibodies protected from SARS-CoV-2 replication in monkeys and mice. Three of 46 monkeys infused with enhancing antibodies had higher lung inflammation scores compared to controls. One monkey had alveolar edema and elevated bronchoalveolar lavage inflammatory cytokines. Thus, while in vitro antibody-enhanced infection does not necessarily herald enhanced infection in vivo, increased lung inflammation can rarely occur in SARS-CoV-2 antibody-infused macaques

Finishing the euchromatic sequence of the human genome

The sequence of the human genome encodes the genetic instructions for human physiology, as well as rich information about human evolution. In 2001, the International Human Genome Sequencing Consortium reported a draft sequence of the euchromatic portion of the human genome. Since then, the international collaboration has worked to convert this draft into a genome sequence with high accuracy and nearly complete coverage. Here, we report the result of this finishing process. The current genome sequence (Build 35) contains 2.85 billion nucleotides interrupted by only 341 gaps. It covers ∼99% of the euchromatic genome and is accurate to an error rate of ∼1 event per 100,000 bases. Many of the remaining euchromatic gaps are associated with segmental duplications and will require focused work with new methods. The near-complete sequence, the first for a vertebrate, greatly improves the precision of biological analyses of the human genome including studies of gene number, birth and death. Notably, the human enome seems to encode only 20,000-25,000 protein-coding genes. The genome sequence reported here should serve as a firm foundation for biomedical research in the decades ahead